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Abstract 

We propose a fast algorithm for computing the expected tranche 
loss in the Gaussian factor model. We test it on portfolios ranging 
in size from 25 ( the size of DJ iTraxx Australia) to 100 (the size of 
DJCDX.NA.HY) with a single factor Gaussian model and show that 
the algorithm gives accurate results. The algorithm proposed here 
is an extension of the algorithm proposed in 4 . The advantage of 
the new algorithm is that it works well for portfolios of smaller size 
for which the normal approximation proposed in [3] in not sufficiently 
accurate. The algorithm is intended as an alternative to the much 
slower Fourier transform based methods El. 



*This work was supported by the Director, Office of Science, Office of Advanced Sci- 
entific Computing Research, of the U.S. Department of Energy under Contract No. DE- 
AC03-76SF00098. 

T E-mail: pokunev@math.lbl.gov 



1 



1 The Gaussian Factor Model 



Let us consider a portfolio of N loans. Let the notional of loan % be equal to 
the fraction fi of the notional of the whole portfolio. This means that if loan 
i defaults and the entire notional of the loan is lost the portfolio loses fraction 
fi or 100/^% of its value. In practice when a loan % defaults a fraction r% of 
its notional will be recovered by the creditors. Thus the actual loss given 
default (LGD) of loan i is 

LGDi — fi(l — n) (1) 

fraction or 

LGDi — 100/j(l - n)% (2) 

of the notional of the entire portfolio. 

We now describe the Gaussian m-factor model of portfolio losses from 
default. The model requires a number of input parameters. For each loan i we 
are give a probability Pi of its default. Also for each i and each k — 1, . . . , m 
we are given a number such that YlT=i w ik < 1- The number w it k is 
the loading factor of the loan % with respect to factor k. Let (pi, ... , <p m and 
(f) 1 , i — 1, . . . , N be independent standard normal random variables. Let <&(x) 
be the cdf of the standard normal distribution. In our model loan i defaults 
if 



W i,k<^k 

k=l 



\ 



i-E«".V <*"'W (3) 

fe=i 



This indeed happens with probability p^. The factors <f>i, . . . , (j) m are usually 
interpreted as the state of the global economy, the state of the regional econ- 
omy, the state of a particular industry and so on. Thus they are the factors 
that affect the default behavior of all or at least a large group of loans in the 
portfolio. The factors 1 , . . . , <p N are interpreted as the idiosyncratic risks of 
the loans in the portfolio. 
Let Li be defined by 

li I{loan i defaulted} (4) 

We define the random loss caused by the default of loan i as 

U = f t (l - r t )I t , (5) 
where is the recovery rate of loan i. The total loss of the portfolio is 

L = J2^ (6) 
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An important property of the Gaussian factor model is that tthe Lj's are 
not independent of each other. Their mutual dependence is induced by the 
dependence of each Li on the common factors (pi, ... , <p m . Historical data 
supports the conclusion that losses due to defaults on different loans are 
correlated with each other. Historical data can also be used to calibrate the 
loadings w^. 

2 Conditional Portfolio Loss L 

When the values of the factors <pi, . . . , m are fixed, the probability of the 
default of loan i becomes 

p* = $-1 ( Pi-X^M (7) 

The random losses Li become conditionally independent Bernoulli vari- 
ables with the mean given by 

E cmd (L i ) = f i (l-r i )p i (8) 

and the variance given by 

VAR cond {Li) = ff (1 - n) V(l - P l ) (9) 

By the Central Limit Theorem the conditional distribution of the portfolio 
loss L, given the values of the factors fa, ... , 4> m , can be approximated by 
the normal distribution with the mean 

E C ond{L) = ^ E C ond{Li) (10) 

i 

and the variance 

VAR cond {L) = ]T VAR cond (Li) (11) 

i 

In j3] it was shown that for portfolios of 125 names this approximation 
leads to accurate results. 

If the size of the portfolio is smaller than 125, for example 30 ( the size 
of DJ iTraxx ex Japan) or 50 (the size of DJ iTraxx CJ), then the Central 
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Limit Theorem no longer provides a sufficiently accurate approximation to 
the conditional distribution of the portfolio loss L. An accurate represen- 
tation of the conditional distribution of the portfolio loss L is given by its 
Hermite series expansion. For historical reasons this expansion is also known 
as the Charlier series expansion CO- 



3 The Hermite Expansion of the Conditional 
Distribution of the Portfolio Loss L 

Let F(x) be the c.d.f. of the conditional distribution of the portfolio loss L. 
So that 

P(L <x) = F(x) (12) 

For each fixed value of the factors (pi, ... , <p m we define the normalized con- 
ditional loss L by 

T = L - E ^ L 1_ (13) 



IVAR cond (L) 

Let F(x) be the c.d.f. of the distribution of the normalized conditional 
portfolio loss L. So that 

P(L <x) = F(x) (14) 
We define the Hermite polynomial H n (x) of degree n by 

2 d n - 2 

#„(*) = (-l)»eV_ e=r (15) 



Let c n be defined by 



Then we have 



in 



H n (x)dF(x) (16) 

-oo 



F(x) = Y: T CiHi{t)^=dt (17) 

oo \/Z7l 



i=0 

The series above converges in the sense of distributions (generalized func- 
tions) [3]. A good reference on the theory of distributions (generalized func- 
tions) is j^J. Let us pick a finite N. Then we have 

N ■ — 

F(x)nJ2ci f Hi(t)^Ldt (18) 

-oo V Z7T 



i=0 



As before the approximation is in the sense of generalized functions. Equation 
(|18|) implies that the distribution of the normalized conditional portfolio loss 
L can be approximated by a distribution with the density 



N 

e 2 



p{x) = y £c i H i (x)-= (19) 

t=0 V27T 

The function p(x) is not necessarily nonnegative and therefore may not be 
a probability density in the strict sense. However, as is explained in [Hj, this 
does not affect the validity of our final result ([24)1 . Therefore we may treat 
p(x) as a real probability density. 

The distribution of the unnormalized loss L can be approximated by a 
distribution with density 

( x ~ E cond( L ) 



p{ X ) = £ , Ci Hi * - EcondiL) 6 - ^ (20) 

i= o JVAR cond (L) \JVAR cond (L) V2n 



The joint distribution of the factors 0i, . . . , <ft m and the portfolio loss L 
can be approximated by a distribution with density 

m 

PjoinMl, ...,<f> m ,L)= p(L) Y[ p G ,Q,l(<frk), (21) 

fe=l 

where Pg.o,i{x) stands for the Gaussian density with mean and variance 1. 

Observe that the coefficient c n depends only on the moments of the dis- 
tribution F(x). Since Lj's are independent Bernoulli random variables these 
moments are known analytically. Thus in the case under consideration all 
the c n 's are known analytically. 

If in equation (|20p we set iV = 1 we obtaine the standard approximation 
by the normal density proposed in jj]. Thus the algorithm proposed here 
is a generalization of the algorithm in 0j. We show later that it gives good 
numerical results even when the portfolio size is too small for the normal 
approximation to be accurate. 
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4 Expected Loss of a Tranche of Loan Port- 
folio 



Let < a < b < 1. We define a tranche loss profile T/ a b (x) by 

mr , . minib — a,max(x — a,0)) . . 

rio,b(x) = — ^ ^ ^ (22) 

Number a is called the attachment point of a tranche, while b is called the 
detachment point of a tranche. The expected loss of a tranche is then 

TLoss(a, b) = J Tl atb (L)p joint ((pi, ...,4> m , . . . 4> m L (23) 

This can be rewritten as a double integral 

TLoss(a, b)= [ f Tl a , b (L)p(L)dL f[ p G , o ,i(0fc)^i ■ • • <t>m (24) 

fc=i 

The inside integral with respect to L can be done analytically for fixed values 
of the factors (f>±, . . . , <p m . The outside integral has to be computed numer- 
ically. However, since it is an integral of a bounded smooth function with 
respect to m-dimensional Gaussian density, it is one of the simpler integrals 
to compute numerically. 

5 Numerical Example 

In this section we apply the proposed algorithm to the single factor Gaussian 
model of a portfolio with n names. We take n to be 25 (size of DJ iTraxx 
Australia), 30 (size of DJ iTraxx ex Japan), 50 (size of DJ iTraxx CJ) and 
100 (size of DJCDX.NA.HY). We choose a single factor model because it is 
the one most frequently used in practice. For each n we compute the loss 
of the equity tranche with the attachment point a = or a = 0% and the 
detachment point 3%. The parameters of the porfolio are 

/. = i 

n 

Pi = 0.015 « 
n — 1 
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r, = OJS- -"*- 1 - 



n — 1 



Wjl = 0.5-H<izi), (25) 
n — 1 

where z = 1, . . . , n. Finally, we choose N = 5 in (|18|). 

In Figure Q we compare the expected loss computed using 10 6 Monte 
Carlo samples with the expected loss computed using formula (J24j) . 1 The 
agreement between the two is good. 



Figure V. Equity Tranche Loss in the Gaussian Single Factor Model 
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6 Conclusions. 

To obtain the results in Figure ^ we only needed to perform a single one 
dimensional numerical integration for each tranche. This is an improvement 

1 The author has the code implementing the algorithm described here in MATLAB, 
VBA for Excel and C. 
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over the Fourier transform based methods |5j which requires computing a 
large number of Fourier transforms for each tranche. Each individual Fourier 
transform is as computationally expensive as (|24j) . 

The expansion (|T5j) is accurate even when the portfolio size is too small 
for the normal approximation of 4J to be precise. Thus we developed an 
algorithm which is as fast as the algorithm proposed in 4J but allows us to 
obtain higher precision for a portfolio of a given size by including more terms 
in ((TBI). 
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